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РЕЗЮМЕ 


В публикации приводятся результаты ГИС-моделирования с помощью программы MaxEnt 3.3.3k области 
распространения средней ящерицы, Lacerta media Lantz et Cyrén, 1920, в прошлом, настоящем и будущем. 
При анализе были использованы географические координаты 26 мест находок вида на территории Дагестана. 
Полученные модели прошлого распространения (последний ледниковый максимум и средний голоцен) в 
северо-восточной части Кавказа демонстрируют стабильное расположение потенциально пригодной терри- 
тории в границах известного современного ареала. Показано отсутствие заметных изменений границ ареала 
вида при потенциальном потеплении климата в будущем (2050 и 2070 гг.). Мы можем сделать вывод об 
отсутствии угроз для средней ящерицы со стороны глобального потепления, кроме того, антропогенное воз- 
действие на дагестанские популяции минимально. Были выявлены наиболее перспективные территории для 
поиска новых популяций: выше и ниже по течению рек Самур, Kypax, Андийское Койсу и Аварское Койсу. 
Кроме того, к ним можно отнести ущелья в среднем течении рек Гамриозень, Халагорк, Уллучай и Чирагчай. 
Отметим, что эта территория слабо изучена герпетологами. 
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ABSTRACT 


The publication presents the results of GIS modelling using the MaxEnt 3.3.3k program for the distribution of the 
Lacerta media Lantz et Cyrén, 1920 in the past, present time and future. During the analysis, the geographical co- 
ordinates of 26 places of finds of the species in the territory of Dagestan were used. The obtained models of past dis- 
tribution (the Last Glacial Maximum and Middle Holocene) in the northeastern part of the Caucasus demonstrate 
a stable location of a potentially suitable territory within the boundaries of the known modern range. It was shown 
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no noticeable changes in the range of the species with potential climate warming in the future (2050 and 2070). We 
can conclude that there are no threats to the species from global warming, in addition, the anthropogenic impact 
on Dagestan populations is minimal. The most promising areas were identified for the search for new populations: 
upstream and downstream the Samur, Kurakh, Andy-Koysu and Avar Koysu rivers. In addition, they include the 
gorges in the middle reaches of the Gamriozen, Halagork, Uluchay and Chiragchay rivers. This territory is poorly 


studied by herpetologists. 


Key words: Lacerta media, area, Dagestan, Caucasus, modeling, MaxEnt. 


BBE/IEHHE 


Средняя ящерица, Гасейа media Lantz et 
Cyrén, 1920 (Рис. 1) — наиболее уязвимый и уз- 
коареальный вид зеленых ящериц рода Lacerta в 
фауне России, занесенный в Красную книгу РФ 
(Даревский [Darevsky] 2001). История изучения 
её распространения началась только в первой по- 
ловине ХХ века после опубликования описания 
Ланца и Сирена. На протяжении длительного 
периода, вплоть до настоящего времени, L. media 
нередко ошибочно определяли как симпатрич- 
ных с ней прыткую, L. agilis Linnaeus, 1758, и по- 
лосатую, L. strigata Eichwald, 1831, ящериц. Это 
значительно затрудняет анализ её распростране- 
ния. Тем не менее за последнее время из печати 
вышел ряд публикаций, посвященных изучению 
вида на северной границе ареала в пределах 
российского Кавказа — Дагестана (Мазанаева и 
Орлова [Mazanaeva and Orlova] 2009; Мазанаева 
и Аскендеров [Mazanaeva and Askenderov] 2016) 
и Краснодарского края (Островских и Плотни- 
ков [Ostrovskikh and Plotnikov] 2009; Туниев и 
Островских [Tuniyev and Ostrovskikh] 2017). 
Необходимо отметить, что эта территория 3a- 
частую не указывается при описании ее ареала 
(Ahmadzadeh et al. 2013b; Smíd et al. 2014). 

Проведение современных исследований ape- 
алов животных уже трудно представить без при- 
менения геоинформационных систем (ГИС). 
Особую популярность приобретает анализ видо- 
вой ниши (Elith et al. 2006; Elith and Leathwick 
2009). Теоретической основой таких работ служит 
гипотеза Гриннелла (Grinnell 1917), что каждый 
вид занимает свою уникальную экологическую 
нишу. Отметим, что параллельно с Дж. Гриннел- 
лом у истоков теории экологической ниши стоял 
и Ч. Элтон; первый из них сближал свой термин с 
местообитанием, поэтому часто к его концепции 
применяют термин «пространственная ниша», 


тогда как второй — с местом вида в пищевой цепи — 
«трофическая ниша». В дальнейшем теория ниши 
получила свое развитие в работах Дж. Хатчинсона 
(Pocheville 2015). Его представление o нише как 
о многомерном объекте, который можно пред- 
ставить в гиперпространстве осей экологических 
параметров, позволило применить для ее ис- 
следования математический анализ, в том числе 
и моделирование (см. подробнее — Лисовский и 
Оболенская [Lissovsky and Obolenskaya] 2014). 

Исследователями была отмечена несомненная 
значимость ГИС-моделирования для изучения 
эктотермных животных, таких как рептилии, с 
ограниченной климатической толерантностью и 
более высокой зависимостью от климатических 
условий (Buckley et al. 2012). Кроме того, по срав- 
нению с птицами и млекопитающими, рептилии 
(и в частности — лацертидные ящерицы) характе- 
ризуются относительно малой мобильностью, что 
также влияет на результативность модели. Ранее 
среди видов Lacertini отечественной герпетофа- 
уны были проанализированы ареалы скальных 
ящериц рода Darevskia с использованием ГИС 
(Доронин [Doronin] 2012, 2015). 

Описание и картирование ареала средней яще- 
рицы на Кавказе было затронуто рядом авторов. 
Представленный в их работах картографический 
материал (процитированы лишь некоторые рабо- 
ты, в основном упоминающиеся в других разделах 
этой публикации) можно разделить по способу 
отображения информации на четыре группы: 
1) нанесение на физическую или административ- 
ную карту известных точек находок ящериц (Алек- 
перов [Alekperov] 1978; Roitberg et al. 2000; Smíd 
et al. 2014; Мазанаева и Аскендеров [Mazanaeva 
and Askenderov] 2016); 2) нанесение предполага- 
емых границ ареалов без конкретизации находок 
(Банников и др. 1971; (Даревский [Darevsky] 
2001; Островских и Плотников [Ostrovskikh and 
Plotnikov] 2009); 3) объединение первого и второ- 
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Рис. 1. Средняя ящерица, Lacerta media Lantz et Cyrén, 1920 (A – половозрелый самец, В — сеголетка), окрестности с. Рутул, 
Рутульский р-н, Дагестан, Россия, 28.07.2009. 

Fig. 1. Three-lined Lizard, Lacerta media Lantz et Cyrén, 1920 (A — ad. male, В — juvenile), vicinities of Rutul, Rutulsky District, Dages- 
tan, Russia, 28.07.2009. 
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ro подхода (Банников и др. [Bannikov et al.] 1977); 
Á) выделение квадратов, построенных Hà основе 
равноугольной графической проекции Меркатора 
(системы ОТМ, 50х50 км), в пределах которых 
известна хотя бы одна точка находки (Sindaco and 
Jereméenko 2008). 

Для Г. media нам известна только одна работа, 
где была использована ГИС-программа — DIVA- 
GIS (Ahmadzadeh et al. 2013а). Согласно данным 
CD. Ахмадзаде с соавторами средняя ящерица име- 
ет наибольшую потенциальную область распро- 
странения среди видов «Lacerta trilineata group»!: 
начиная от Загросских гор в западном Иране, до 
центральных районов Греции, что значительно 
превосходит реальные границы ареала вида. Это 
соответствует области распространения трех- 
линейчатой ящерицы, L. trilineata Bedriaga, 1886, 
согласно воззрениям на объем вида в середине 
ХХ века, т.е. когда L. media рассматривалась как ee 
подвид. Однако в этой публикации для всего рос- 
сийского Кавказа, где проходит северная граница 
ареала номинативного подвида, были учтены 
только три точки, что не могло не повлиять на до- 
стоверность полученных данных. Неудивительно, 
что на полученных картах северный макросклон 
Кавказа полностью выпал из потенциального рас- 
пространения средней ящерицы. При этом Кав- 
каз — перспективная территория для филогенети- 
ческих исследований зеленых ящериц, поскольку 
здесь предполагается наличие целого ряда рефу- 
гиумов различного возраста, в том числе и северо- 
восточного Кавказа — семиаридные котловины 
внутригорных районов Дагестана (Мазанаева и 
Туниев [Mazanaeva and Tuniyev] 2011). 

Не так давно были разработаны новые под- 
ходы поиска рефугиумов, связанные с ГИС-мо- 
делированием ареалов (Скоринов и Литвинчук 
[Skorinov and Litvinchuk] 2016). Основным по- 
ложением, на котором они базируются, выступает 
предположение, что виды характеризуются ста- 
бильными нишами, не меняющимися, по крайней 
мере, в течение последних нескольких десятков 
тысяч лет. Как следствие, изучив современные 
климатические предпочтения вида, можно экс- 
траполировать эти данные на прогнозируемые 
условия прошлого, что становится дополнитель- 
ным источником информации (в совокупности 
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с палеонтологическими и филогенетическими 
данными) для изучения истории формирования 
ареалов. Это же применимо и к моделированию 
их области распространения в будущем, что ис- 
пользуется (в том числе и в герпетологических 
работах) для выявления динамики потенциаль- 
ного изменения ареала и оценки предполагаемого 
ущерба изученному виду от последствий глобаль- 
ного потепления (Литвинчук [Litvinchuk] 2017). 

Исходя из вышесказанного, целью данной 
работы было создать модель прошлого, совре- 
менного и будущего ареала средней ящерицы на 
территории Дагестана, дать характеристику рас- 
положения потенциально пригодных для нее тер- 
риторий, что имеет высокую природоохранную 
значимость. 


МАТЕРИАЛЫ И МЕТОДЫ 


Материалами для данного исследования по- 
служили коллекции Зоологического института 
РАН (Санкт-Петербург, ЗИН РАН), Зоологи- 
ческого музея Московского государственного 
университета (Москва, ЗМ МГУ), Харьковского 
национального университета им. В.И. Каразина 
(Харьков, МП ХНУ), собственные полевые ис- 
следования 1980—2018 гг., а также литературные 
данные. 

Всего наш список находок £L. media на террито- 
рии Дагестана включил 26 локалитетов (Табл. 1). 
Определение географических (десятичных) 
координат в полевых условиях производили с 
помощью GPS навигатора (Garmin); при работе 
C музейными коллекциями и литературными ис- 
точниками использовали интерактивные карты 
(www.wikimapia.org и www.maps.google.ru) 

При построении моделей использовали 19 
биоклиматических параметров, отражающих 
данные по температуре и осадкам территории 
(Busby 1991), которые позволяют провести ин- 
терполяцию наблюдаемых данных с 1950 по 
2000 rr. (Current) (Hijmans её al. 2005). Эти по- 
казатели имеют наибольшее значение для рас- 
пределения наземных биологических объектов 
(Hijmans et al. 2005). Мы не использовали данные 
по высоте местности, так как от них в значитель- 


По нашему мнению здесь необходимо использовать термин «надвидовой комплекс» (superspecies complex) и 


обозначение Lacerta (trilineata). 
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ной степени линейно зависимы температурные 
параметры (Harris et al. 2013). Кроме того, при 
моделировании ареала средней ящерицы с дан- 
ными по высотам и без них полученные модели 
распространения оказались идентичными. 

Этисведения были получены избазы Worldclim 
(http://www.worldclime.org) в разрешении 30 arc- 
секунд и 10 arc-minute. В этой же базе были взяты 
слои, описывающие климатические условия во 
время последнего ледникового (гляциального) 
максимума (около 22 тыс. лет назад) (Last Glacial 
Maximum) и среднего голоцена (около 6 тыс. лет 
назад) (Mid-Holocene). 

С помощью программы ENM Tools 1.4.3 (War- 
ren et al. 2010) были удалены локалитеты, близко 
расположенные друг к другу и находящиеся в 
одном и том же квадрате (1X1 км) и проведен 
корреляционный анализ между слоями. При ана- 
лизе хорологии ящериц использовали программу 
MaxEnt (Maximum Entropy Species Distribution 
Modelling, версию  3.33k) (wwwes.princeton. 
edu/-schapire/maxent/), в которой экстраполяция 
(моделирование) географического распростра- 
нения биологических видов проводится методом 
максимальной энтропии. Метод максимальной 
энтропии применяется для поиска свойств рас- 
пределения факторов среды в наборе точек про- 
странства (ячеек изучаемого растра), в которых 
вид обитает: по исходной гипотезе в таком наборе 
свойства распределений факторов среды, вли- 
яющих на распространение вида, должны быть 
максимально близки свойствам распределений 
этих факторов во всем изучаемом географическом 
пространстве (Elith et al. 2011). 

C помощью данной программы созданы карты 
наиболее вероятных областей распространения 
ящериц и определен вклад каждого фактора B по- 
строение модели. Программа МахЕп позволяет 
определить степень влияния параметров среды 
(в %) на границы области распространения из- 
учаемого таксона (Phillips and Dudik 2008). Для 
тестирования полученных моделей использовали 
25% точек. На начальном этапе при построении 
первичной модели производилась оценка значи- 
мости биоклиматических параметров, а впослед- 
ствии — исключение малозначимых и коррелиру- 
ющих друг с другом. В результате этих действий 
оставлялась группа из 3—4 наиболее значимых 
параметров, на основе которых строились окон- 
чательные модели. Оценка достоверности модели 


проводилась с использованием встроенной функ- 
ции построения кривой AUC (area under receiver 
operating characteristic (ROC) curve), отражающей 
чувствительность вида к параметрам и ero спец- 
ифичность. Значение индекса AUC выше 0.75 и 
близкое к единице означает высокую степень до- 
стоверности построенной модели (Elith 2002). 
Для оценки того, как климатические измене- 
ния повлияют на пространственное распределение 
L. media в будущем (2050 и 2070 rr.), также исполь- 
зовали растровые слои из базы WorldClim, рассчи- 
танные с использованием хорошо зарекомендо- 
вавшей себя (Brandefelt et al. 2013) климатической 
модели CCSMA (Gent et al. 2011) для четырёх pe- 
презентативных траекторий концентраций (RCP), 
разработанных Межправительственной группой 
по изменению климата (IPCC). КСР являются 
сценариями климатических состояний, характе- 
ризующими величину антропогенно-обусловлен- 
ного радиационного воздействия, достигаемого к 
2100 г. по сравнению c 1750 r. (2.6, 4.5, 6.0 и 8.5 Вт/ 
м?). Согласно этим сценариям вероятные оценки 
увеличения глобальной температуры к концу 
2100 г. составят: 0.2—1.8 °С (ВСР 2.6); 1.0-2.6 °С 
(ВСР 4.5); 1.3-3.2 °С (RCP 6.0); 2.6—4.8 °С (RCP 
8.5). Визуализацию ГИС-карт провели с помощью 
программы DIVA-GIS 7.5.0 (www.diva-gis.org) 
(Scheldeman and van Zonneveld 2010). 


РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 


В нашем исследовании при моделировании 
прошлой, настоящей и будущей области распро- 
странения вида были получены хорошие резуль- 
таты производительности модели потенциального 
распределения (индекс AUC > 0.998). Полученные 
карты по этим данным надежно характеризуют 
особенности распространения изученного вида и 
позволяют нам в совокупности с новым материа- 
лом дать свои замечания и дополнения. 

Наиболее значимый вклад в построение моде- 
ли (> 10 %) современного ареала вида внесли (в 
порядке убывания) количество осадков наиболее 
жаркого сезона (bio 18), средняя температура 
наиболее холодного сезона (bio 11), количество 
осадков за сезон (bio 15) и средняя температура 
наиболее сухого сезона (bio 9) (Табл. 2). Высо- 
кий вклад параметров bio 9 и 11 характеризуют 
L. media как обитателя территории с выраженным 
континентальным климатом; с недостаточным 


ATO 


увлажнением (400—600 MM B год) можно связать 
доминирование bio 15 и 18. 

В исследованиях Ф. Ахмадзаде с соавторами 
(Ahmadzadeh et al. 2013a) наибольший вклад в 
построение глобальной модели ареала средней 
ящерицы внесли сезонная температура (bio 4), 
годовое количество осадков (bio 12), количество 
осадков наиболее влажного сезона (bio 16) и 
количество осадков наиболее холодного сезона 
(bio 19). Эти показатели полностью не соответ- 
ствуют нашим данным, что, по-видимому, может 
свидетельствовать о значительной региональной 
специфике условий обитания вида, особенно в 
условиях горного ландшафта. 

Наиболее значимый вклад в построение мо- 
дели среднеголоценового ареала вида внесли ко- 
личество осадков за сезон (bio 15), среднегодовая 
температура (bio 1), минимальная температура 
наиболее холодного месяца (bio 6), количество 
осадков наиболее сухого месяца (bio 14). Они во 
многом соответствуют показателям модели ареа- 
ла в гляциальный максимум: также преобладают 
bio 1, 6, 15; отличие заключается в доминировании 
bio 9 и количества осадков наиболее сухого сезона 
(bio 17). 

При моделировании вероятного ареала с уче- 
том климатических трендов к 2050 и 2070 rr. на- 
бор ключевых параметров оказался практически 
идентичным при всех сценариях потепления — bio 
9, 11, 14, 15. Эти же четыре параметра выходят 
на первый план при моделировании ареала про- 
шлого и настоящего, что говорит об их стабиль- 
ном вкладе в формирование «пространственной 
ниши» средней ящерицы. 

В публикациях по моделированию ареалов 
различных групп организмов неоднократно была 
отмечена связь широты экологической амплитуды 
видов с качеством моделирования и структурой 
вкладов используемых параметров: чем вид более 
стенотопен, тем меньше набор параметров окружа- 
ющей среды, доминирующий по вкладу в постро- 
ение ГИС-модели ареала (см. подробнее — Дудов 
[Dudov] 2016). В нашем случае превалирование 
четырех параметров и наличие среди них одного 
ведущего говорят o стенотопности L. media. 

На ГИС-модели современного распростране- 
ния вида (Рис. 2) наиболее подходящая с клима- 
тической позиции область выявлена в ландшафт- 
ной зоне передовых хребтов Дагестана, отчетливо 
формирующих северную и восточную границу 
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ядра ареала, внутригорной и высокогорной зо- 
нах — южная и западная границы. Мы можем 
ожидать новых находок средней ящерицы выше 
и ниже по течению рек Самур, Kypax, Андийское 
Койсу и Аварское Койсу. Кроме того, перспектив- 
ными являются ущелья в среднем течении рек 
Гамриозень, Халагорк, Уллучай и Чирагчай. Отме- 
тим, что эта территория слабо изучена герпетоло- 
гами. Интересно сопоставить ГИС-ранжирование 
территории по пригодности с данными полевых 
исследований: в районе с. Ботлих плотность на- 
селения вида составляет 7.5—10 экз./га, при этом 
территория отнесена к зоне средней (74-56%) и 
низкой (56-37%) вероятности обитания, в то вре- 
мя как в районе сел Ахты, Рутул и Тидиб средняя 
ящерица более многочисленна (12—62 экз./га), а 
территория вошла в ядро ареала с максимальной 
вероятностью (100-74%). 

Согласно ранее опубликованным данным 
(Мазанаева и Аскендеров [Mazanaeva and Asken- 
derov] 2016) дагестанские популяции изолирова- 
ны от закавказской части ареала вида. Это делает 
данную географическую популяцию удобным 
объектом для ГИС-моделирования. В то же время 
на картах в определителях А.Г. Банникова с со- 
авторами (Банников [Bannikov] 1971, 1977) и в 
очерке А.Л. Агасяна с соавторами (Agasyan et al. 
2009) дизъюнкции (по крайней мере, в долине 
р. Самур) отсутствуют. Полученные полевые дан- 
ные и ГИС-карты говорят в пользу наличия изо- 
ляции дагестанских популяций: на картах четко 
прослеживается южная граница благоприятной 
области для этой ящерицы в регионе. Это COOT- 
ветствует изоляции Северо-Восточного и Куро- 
Араксинского рефугиумов средиземноморской 
герпетофауны на Кавказе (Tuniyev 1995). 

В понимании Н.П. Наумова (Наумов [Nau- 
mov] 1963) географическая популяция «населяет 
территорию с географически однородными ус- 
ловиями существования — зону сходного благо- 
приятствования» (с. 307). Это соответствует BHI- 
явлению на полученных ГИС-картах отдельных 
ядер ареала – территорий с максимально под- 
ходящими условиями и вероятностью обнаруже- 
ния на ней изучаемого объекта. Географическая 
популяция, в свою очередь, распадается на более 
мелкие пространственные группировки — эколо- 
гические популяции, состоящие из обитателей 
одного местообитания (биотопа). В нашем случае 
под это понятие подходят популяции L. media от- 
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Рис. 2. Модель современного потенциального распространения средней ящерицы, Lacerta media, созданная с помощью IIporpaM- 
мы MaxEnt 3.3.3k. Использованы только точки находок с территории Дагестана (показаны кругами черного цвета). Внизу — уве- 
личенный фрагмент изученной территории. Здесь и далее: цветная заливка обозначает области вероятного распространения с 
высокой (100—74%; красный цвет), средней (74-56%; оранжевый цвет), низкой (56-37%; желтый цвет) и крайне низкой (37-19%; 
зеленый цвет) вероятностью встречи. 

Fig. 2. Potential distribution modeling of Lacerta media obtained with MaxEnt 3.3.3k. Used only localities in Dagestan (marked by black 


circles). Below is an enlarged fragment of the studied territory. Here and below: colors in the map designate different suitability values: 
high (100—74?6; red), medium (74—56?6; orange), low (56-37%; yellow), and extremely low (37-19%; green). 
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дельных семиаридных котловин крупных горных 
рек (Рис. 3). 

В предшествующей публикации (Мазанаева 
и Аскендеров [Mazanaeva and Askenderov] 2016) 
ставилось под вопрос обнаружение средней яще- 
рицы в предгорье Главного Кавказского хребта 
между г. Варташеном (= г. Огуз) и с. Халкалы 
Азербайджана на южном склоне Главного Кав- 
казского хребта (Алекперов [Alekperov] 1978). 
До публикации Алекперова эта территория 
была включена в область распространения 
ящерицы на карте в монографии Банникова 
с соавторами ([Bannikov et al.] 1977). Данные 
ГИС-моделирования говорят о наличии здесь 
подходящих климатических условий для вида. 
Область потенциального распространения вида 
(но с меньшей вероятностью обнаружения) 
простирается и на сопредельную с Дагестаном 
территорию Чечни (Веденский, Ножай-Юртов- 
ский, Шатойский р-ны), где имеются аридные 
котловины востока Северо-Юрской депрессии. 
Они также были отнесены Б.С. Туниевым к 
ксерофитным рефугиумам, отличающимся друг 
от друга набором видов геретофауны (Tuniyev 
1995). Карта потенциальных ареалов, получен- 
ная нами с использованием программы MaxEnt, 
в общих чертах верно характеризует особенности 
распространения вида, однако в ряде случаев 
она содержит явные искажения, поэтому ее про- 
гностическая роль не должна переоцениваться. 
Так, в благоприятную область попали не при- 
годные для вида территории — высокогорье на 
севере Азербайджана и юго-западе Дагестана 
(Тляратинский р-н). Как было отмечено нами 
при работе с моделированием ареалов видов 
герпетофауны Крыма (Кукушкин и Доронин 
[Kukushkin and Doronin] 2015), по-видимому, 
сетка базы WordClim c ячейками 1х1 км является 
слишком грубой для сравнительно небольших, 
но сложных в отношении рельефа и климатиче- 
ских условий таких горных регионов, как Кавказ 
и Крым, не отражает в полной мере разнообразие 
мезоклиматов и не учитывает микроклимати- 
ческие особенности территории. Помимо того, 
что программа в силу своих характеристик не 
способна учесть все многообразие факторов, 


обусловливающих распространение животного, 
на результаты оказывает влияние неточность 
климатических баз, неравномерность изученно- 
сти ареала и его сохранность, а также, вероятно, 
история заселения видом региона (т.е. несоответ- 
ствие фактического и потенциального ареалов). 

Еще одной территорией, отмеченной на всех 
полученных картах как благоприятной для оби- 
тания средней ящерицы, является регион Кав- 
казских Минеральных Вод, где средняя ящерица 
достоверно неизвестна, но проходят границы 
ареалов многих средиземноморских видов, либо 
имеются их изолированные участки ареала. Так, 
по данным А.Д. Михеева (Михеев [Mikheev] 
2005) подавляющее большинство реликтов во 
флоре этой территории средиземноморского про- 
исхождения. Из списка фауны КМВ наиболее 
наглядный пример — оливковый полоз, Platyceps 
najadum (Eichwald, 1831) (вероятнее всего — вы- 
мерший здесь вид), ставший вместе со средней 
ящерицей индикатором наличия рефугиумов 
средиземноморских видов герпетофауны Кавказа 
(Tuniyev 1995). 

Обе полученных ГИС-модели прошлого рас- 
пространения в северо-восточной части Кавказа 
(Рис. 4) демонстрируют стабильное расположение 
потенциально пригодной территории (100—56% 
вероятности обнаружения) в границах известного 
современного ареала. Это еще раз подтверждает 
положение о наличии на данной территории ре- 
фугиума. Как и следовало ожидать для этого ксе- 
рофильного вида, площадь пригодной территории 
во время последнего ледникового (гляциального) 
максимума заметно уступает таковой на среднего- 
лоценовой модели. 

Наши интерпретации полученных карт пале- 
оареалов не могут быть подтверждены или опро- 
вергнуты фактическими палеонтологическими 
находками в силу их отсутствия на территории 
Дагестана. Единственное указание на подобную 
находку в пределах Кавказа — обнаружение ниж- 
ней челюсти средней ящерицы в мустьерских 
слоях (верхний плейстоцен) пещеры Швалиети 
Грузии (Векуа [Vekua] 1991). При этом указан- 
ный материал требует более детального изучения 
и подтверждения видового определения”. По 


Несомненный интерес представляют и другие находки Lacertidae из плейстоцен-голоценовых местонахож- 
дений Грузии с отсутствующим родовым и (или) видовым определением (Буачидзе [Buachidze] 1973; Алифанов 


[Alifanov] 2012; Blain et al. 2014). 
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Рис. 3. Места обитания средней ящерицы, Lacerta media: A — окрестности c. Ахты, ущелье p. Ахтычай, 4.05.2013; В — ущелье 
р. Самур, 2.05.2012. 


Fig. 3. Habitat of Lacerta media: А — the vicinity of Akhty, the gorge of the Akhtychay River, 4.05.2013; B — the gorge of the Samur River, 
2.05.2012. 
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Рис. 4. Модели потенциального распространения средней ящерицы, Lacerta media во время последнего ледникового (гляциально- 
го) максимума (А) и среднего голоцена (В), созданные с помощью программы MaxEnt 3.3.3k. 


Fig. 4. Potential distribution modeling of Lacerta media in the Last Glacial Maximum (A) and Mid-Holocene (B) obtained with MaxEnt 
3.3.3k. 
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Рис. 5. Модели потенциального распространения средней ящерицы, Lacerta media к 2050 г., созданные с помощью программы 
MaxEnt 3.3.3k: A — RCP 2.6; В - RCP 4.5; С - ВСР 6.0; D - RCP 8.5. 


Fig. 5. Potential distribution modeling of Lacerta media by 2050 obtained with MaxEnt 3.3.3k.: A – RCP 2.6; B - RCP 4.5; С - RCP 6.0; 
D - RCP 8.5. 
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Рис. 6. Модели потенциального распространения средней ящерицы, Lacerta media к 2070 r., созданные с помощью программы 
MaxEnt 3.3.3k: A – RCP 2.6; B – RCP 4.5; C – ВСР 6.0; D - ВСР 8.5. 


Fig. 6. Potential distribution modeling of Lacerta media by 2070 obtained with MaxEnt 3.3.3k.: A – RCP 2.6; B - RCP 4.5; C - RCP 6.0; 
D - RCP 8.5. 
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той же причине мы не можем привлечь и сведе- 
ния по филогении дагестанских популяций и 
датировкам их «молекулярных часов». Но наши 
данные не противоречат датировкам для вида в 
целом (Ahmadzadeh et al. 2013b): отделение клады 
«Lacerta media» в «Lacerta trilineata group» произо- 
шло B плиоцене (около 3.5 млн. лет назад) и было 
спровоцировано горообразованием в Анатолии, а 
дальнейшие климатические колебания сформи- 
ровали современную картину расхождения ли- 
ний. Появление закавказской субклады датируют 
1 млн. лет, при этом расселение средней ящерицы 
в регионе шло, вероятнее всего, из Гирканского 
рефугиума (в указанной статье он назван Эльбур- 
ским). 

Карты будущего распространения средней 
ящерицы (Рис. 5, 6) демонстрируют в целом ста- 
бильность расположения и площади пригодных 
местообитаний в Дагестане. Мы можем сделать 
вывод об отсутствии угроз для вида со стороны 
глобального потепления. Основная причина со- 
кращения численности вида на Черноморском 
побережье Краснодарского края – фрагментация, 
трансформация и уничтожение местообитаний, 
пресс браконьерства (Туниев и Островских 
[Tuniyev and Ostrovskikh] 2017), в то время как 
антропогенное воздействие на дагестанские по- 
пуляции минимально. 

Дальнейшие исследования филогении и фило- 
географии кавказских популяций L. media позво- 
лят оценить полученные с помощью ГИС данные, 
более детально проследить историю расселения 
вида и формирование современных рефугиумов. 
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